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Driving spin excitations by hydrostatic pressure in BiFeOa 
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Optical spectroscopy has been combined with computational and theoretical techniques to show 
how the spin dynamics in the model multiferroic BiFeOa responds to the application of hydrostatic 
pressure and its corresponding series of structural phase transitions from R3c to the Pnma phases. 
As pressure increases, multiple spin excitations associated with non-collinear cycloidal magnetism 
collapse into two excitations, which show jump discontinuities at some of the ensuing crystal phase 
transitions. Effective Hamiltonian approach provides information on the electrical polarization and 
structural changes of the oxygen octahedra through the successive structural phases. The extracted 
parameters are then used in a Ginzburg-Landau model to reproduce the evolution with pressure of 
the spin waves excitations observed at low energy and we demonstrate that the structural phases 
and the magnetic anisotropy drive and control the spin excitations. 

PACS numbers: 77.80.Bh, 75.50.Ee, 75.25.+Z, 78.30.Hv 


Multiferroic insulators have non-collinear magnetic or¬ 
der that drives ferroelectricity, and ferroelectric order 
that controls magnetism. This demonstrated cross¬ 
correlation between electric and magnetic effects shows 
great promise for the development of magnonic devices 
whose goal is to use magnetic excitations as a low en¬ 
ergy substitute of conventional electronics [iHS]. The 
spin-lattice interaction plays a decisive role in mediating 
the combined ferroic properties of multiferroic materials. 
When induced by epitaxial mismatch or chemical sub¬ 
stitutions, strain provides a handle for the complex in¬ 
terplay between magnetic and electronic properties and 
their coupling to structural distortions ms]. At present, 
the mechanisms linking spin excitation to structural de¬ 
formation remain hardly accessible, whereas their fine 
control is highly desirable to build these new technolo¬ 
gies. 

Bismuth ferrite (BiFeOs) plays an important role in 
multiferroics research as it is one of the few materi¬ 
als that has coexisting ferroelectricity and magnetism 
at room temperature [7] with an unusual combination 
of properties such as large above band gap voltages 
[5], photovoltaic effect [3] and conductive domain walls 
m- At ambient pressure, it becomes ferroelectric below 
~ 1100 K, with one of the largest known electrical po¬ 
larization (that is, P = 100 fiC/cm^). Below 640 K, it 
exhibits an antiferromagnetic-like spiral of the cycloidal 


type with wavevector Qq = 27r/64nm |llj . This spiral 
transforms into a canted homogenous antiferromagnetic 
state under epitaxial strain 0 i or chemical doping 
[^. Hydrostatic pressure is known to induce dramatic 
changes in BiFeOa’s ferroelectricity and crystal struc¬ 
ture. The large electric polarization either disappears or 
becomes weak above « 5 GPa, and a total of six struc¬ 
tural phase transitions have been observed up to 50 GPa 

[ung. 

In this Letter, we combine an advanced high-pressure 
technique ideally suited to probe simultaneously spin and 
phonon excitations with Landau-Ginzburg and effective 
Hamiltonian calculations to elucidate the coupling be¬ 
tween spin excitations and structure in the prototypi¬ 
cal multiferroic BiFeOa. We determine how pressure- 
induced structural transitions drive the magnetic order 
from a non-collinear to an homogeneous magnetic state. 

We performed Raman spectroscopy measurements on 
single crystals of bulk BiFeOa under hydrostatic pressure 
up to 12 GPa in a membrane diamond anvil cell. We have 
developped an original optical experimental setup in or¬ 
der to track low energy excitations down to 7 cm“^ un¬ 
der extreme conditions (See Supplemental Material). We 
have thus been able to follow simultaneously the phonon 
modes and the magnetic excitations under a broad range 
of hydrostatic pressure. 

Figure]^ a) shows the phonon modes and the spin exci- 
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FIG. 1: a) Raman spectrum of spin wave excitation modes 
at low energy in BiFeOs under 0.99 GPa and at 300 K. The 
assignment of these spin excitation modes is in agreement 
with Refs. [261 [27] . The first spin wave excitation peak at 
7 cm“^ (denoted by ”?”) has not been assigned (See text). 
Inset: Large energy scale Raman spectrum of phonon modes 
b) Magnetic incommensurate cycloid in BiFeOa and the two 
sets of spin wave excitations (p and ip as in-plane cyclon and 
extra-cyclon modes, respectively. The electrical polarization 
vector P is along [111] and the cycloid propagating along 
[1,0,-lj lies in the (-1,2,-!) plane. 

tations measured at low pressure. The frequency and the 
optical selection rules of the phonons modes are charac¬ 
teristic of the rhombohedral (R3c) phase [14] . The series 
of narrow peaks (with linewidth ^ 1 cm“^) in Fig. 
a) are the fingerprint of the cycloidal spin excitations at 
zero wavevector m- The spin excitations in BiFeOa can 
be decomposed into cyclon (^„) and extra-cyclon {ipn) 
modes, corresponding to oscillations in and out of the cy¬ 
cloid plane, respectively (Fig. 0 b) UMI]. The first spin 
wave excitation peak at 7 cm“^ has not been assigned 
because it can either be attributed to the mode or 
can result from small domains in the sample with weak 
magnetization m- 

As indicated by the change in the phonon modes 
(See Fig 0 a and Supplemental Material), we observe 
four structural transitions occurring at about 3.5, 5.5, 
7.75 and 11 GPa, from a rhombohedral (R3c) to an or¬ 
thorhombic (Pnma) phase through three orthorhombic 
structures (Oi, O 2 , O 3 ), in agreement with previous high 
energy Raman scattering and X-ray studies [T^ . Thanks 
to our observation of the new phonon modes occurring 


below 100 cm all structural transitions can be tracked, 
especially the transitions between the O 2 and O 3 phases. 




FIG. 2: a) Evolution of the energies of selected phonon modes 
in the 40-200 cm“^ range. Several modes appear or disappear 
at the structural transitions, b) Effective Hamiltonian simu¬ 
lations of the structure : the plotted lines corresponds to the 
enthalpies of four different structures as a function of pres¬ 
sure, relative to the enthalpy of the R3c phase. As pressure 
increases we observe the sequence of phase transitions R3c —>■ 
Pna2i —> Pca2i —>■ Pnma. 

To support the structural description of BiFe 03 un¬ 
der pressure, we performed a theoretical study of the 
crystalline structure using the effective Hamiltonian ap¬ 
proach developed in Refs. [m HOj. As shown in Fig. 0 
b), at low pressure the crystalline structure is rhombo¬ 
hedral and belongs to the R3c space group; when the 
pressure increases, we find several transitions towards or¬ 
thorhombic structures with complex oxygen octahedra 
tilts belonging to the space groups Pna2i and Pca2i. 
These complex structures belong to the family of nan¬ 
otwin phases predicted in Ref. |2D| that have energies 
close to those of the R3c and Pnma states. Finally, 
at the highest pressures, theory predicts a transition to 
the orthorhombic structure belonging to the Pnma space 
group. One can thus notice that all but one of the mea¬ 
sured structural phases are reproduced in the calcula¬ 
tions. This missing phase is likely another intermedi¬ 
ate and stable nanotwin configuration that may have a 
slightly higher enthalpy in the present effective Hamil¬ 
tonian calculations. This theoretical study also provides 
information on the electrical polarization and structural 
changes of the oxygen octahedra through the successive 
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FIG. 3: a) Low energy part of some of the Raman spectra showing the behaviour of the magnetic excitations under pressure, 
b) Energy of the spin wave excitations in BiFeOs from 0 to 12 GPa reported as a function of pressure. Lines are fits using our 
theoretical model (See text) and some structural parameters obtained from computations. Vertical dash lines mark the four 
structural transitions. 


structural phases. These extracted parameters will be 
shown to be crucial to reproduce the evolution with pres¬ 
sure of the spin waves excitations observed at low energy. 


Figure a) shows the low-energy part of the Raman 
spectra obtained at different pressures, and Fig.|^b) de¬ 
picts the corresponding spin excitation energies as a func¬ 
tion of the applied pressure. At low pressure, below the 
first structural transition at 3.5 GPa, we observe that 
some spin excitations harden while other soften with a 
clear tendency for the three lowest energy modes to merge 
towards 12 cm“^. Otherwise, the width of these peaks 
remains constant under pressure which indicates good hy- 
drostaticity. Above 3.5 GPa, the crystal structure enters 
the first orthorhombic phase and only two spin excita¬ 
tions are observed, signaling the sudden disappearance 
of the spin cycloid at the first structural transition. The 
presence of two spin excitations in the Oi phase shows 
that in this pressure range, BiFeOa’s magnetic order is 
a simple two-sublattice antiferromagnet. The spin exci¬ 
tations harden as the pressure increases. A jump dis¬ 
continuity is observed at the Oi-to -02 phase transition 
whereas the hardening is continuous through the 02 -to- 
O 3 phase transition. 


In order to describe these results we propose a modified 
model for the coupling between magnetism and ferroelec- 


tricity in BiFeOa [lain] 

H = ^ X! I X 5)^D{u}i- Wi+g)] 

i,S I 

• 5 ,)^ ( 1 ) 

with the sum running over all sites i of the pseudo cubic 
lattice formed by the Fe^'*' ions, with 8 being the unit 
vector linking each site to its six nearest neighbours. The 
vector Si describes the Fe^+ spin at site 1 , while the 
vectors Ui and Ui represent structural distortions around 
this site. More precisely and as schematized in Fig. Ui 
is the local mode describing the electric dipole moment at 
site i (the ionic contribution to the electrical polarization, 
P, is proportional to (wi)), and uii is a pseudovector 
describing oxygen octahedra tilting at site i (the direction 
of cUi is the axis about which the oxygen octahedron of 
site i tilts, while its magnitude is the value of the angle 
associated with such tilt | 22 ) . 

All interaction energies in Eq. 0 are expected to be 
pressure dependent. The exchange interaction J > 0 is 
known to scale inversely proportional to the tenth power 
of the distance between the Fe spins [3S] , and to increase 
linearly with pressure in each structural phase [12]. The 
pressure dependence of the other model parameters is not 
known. The parameter C is the spin-current interaction 
responsible for cycloidal order in BiFeOa, and the param¬ 
eter D leads to spin canting and weak ferromagnetism 
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FIG. 4: a) Calculated average local mode (n) (directly pro¬ 
portional to P) for each phase as a function of pressure. Inset: 
schematic representation of the local mode u, the polarization 
P and the antiferrodistortive vector ujafd describing the tilt¬ 
ing of BiFeOa’s oxygen octahedra in the pseudo-cubic struc¬ 
ture of BiFeOs. The Dzyaloshinskii-Moryia interaction vector 
points along (jJafd- b) Calculated antiferrodistortive vector 
ij^AFD (left scale) and the 8 angle between (u) and ijJafd 
(right scale) as a function of pressure. 


[HIM]- Both C and D arise from the Dzyaloshinskii- 
Moriya interaction, with C,D o: rjsoJ^ with riso being 
the spin-orbit energy splitting of the heaviest lattice ion, 
Bi^+. Finally, the last term of Eq. ([^ models single ion 
anisotropy, with the anisotropy axis pointing along the 
local mode vector with K oc rjgQ as shown in a mi¬ 
croscopic calculation m- 

At ambient pressure and room temperature, BiFeOs 
has the R3c structure with the same local mode, u, at 
all sites i, pointing along [111] and giving rise to P. In 
addition, the oxygen octahedra adopts an antiferrodis¬ 
tortive order, with iVi = in one sublattice and 

oUi+g = —loafdI^ in the other, with ujafd also pointing 
along [111]. This structural configuration combined with 
Eq. (0 gives rise to the magnetic cycloidal order shown 
in Fig. [^b). When the effective anisotropy 

= Ku^ - ^^^afd ( 2 ) 

is positive, the ground state of Eq. 0 is an anharmonic 


cycloid, with spin excitations splitting into waves of even 
(-I-) and odd (—) symmetry with respect to space inver¬ 
sion along the cycloid wavevector Q (See Supplemental 
Material). This splitting is found to be essential to de¬ 
scribe the modes [55H1Z]; for example, at ambient pres¬ 
sure, the modes are split by 5.5 cm’ -1 m- From our 
model calculations, we expect the evolution with pressure 
of the splitted modes to track the pressure dependence 
of J, u, and uiafd (See Supplemental Material). 

Our effective Hamiltonian calculations predict that u 
decreases with increasing pressure and ujafd remains 
nearly constant in each phase (Fig.0, as consistent with 
the known fact that (small) pressure typically reduces 
ferroelectricity [28]. In Fig. [^b), the pressure trend of 
the spin excitations is well reproduced by our simple two 
parameters model. Only one of them depends on pres¬ 
sure (See Supplemental Material). 

With pressures above about 3.5 GPa, a structural 
phase transition has occured, and only two spin exci¬ 
tations with frequencies to be denoted by Whigh and Wiow 
are observed. Equation 0 predicts a transition to a 
(canted) homogeneous antiferromagnet when 


^ JlATeffl _ i^AFD 

8 {CuY J 


(3) 


so we infer that either J decreases or ujafd/ u increases 
during the transition R3c —?► Oi. Both behaviours are 
consistent with our effective Hamiltonian calculations. 
Obviously if P becomes zero right at the Oi phase bound¬ 
ary, a transition towards an homogeneous antiferromag¬ 
netic state would appear associated with a single spin 
waves peak. This scenario is inconsistent with the mea¬ 
surements and the structural calculations. Indeed, u is 
not zero in the nanotwin phases as evidenced in Fig. [^ 
a). We combine the two theoretical models to reproduce 
the spin wave excitations in the orthorhombic phases. 
Equation 0 is able to explain the observed magnon 
data provided that two conditions are satisfied: K < 0 
and 0 yl 0 , where 9 is the angle between (u) and ujafd 
shown in Fig. Only under these two conditions do we 
get Whigh > v 2 wiow with wiow > 0 as observed experi¬ 
mentally in the homogeneous magnetic phases. 

Turning again to our effective Hamiltonian calcula¬ 
tions, we see that 9 becomes nonzero only after the tran¬ 
sition R3c —>■ Pna2i; and that 9 has a jump disconti¬ 
nuity in the transition Pna2i —>■ Pca2i (Fig. 0. We 
do not define 9 for the non-polar phase Pnma. Using 
the effective Hamiltonian calculations of 9 for Pna2i and 
Pca2i to model magnons in the Oi and O 2 -O 3 phases, 
respectively, we are able to reproduce the pressure de¬ 
pendence of the two spin wave excitations in the Oi, O 2 
and O 3 (Fig. 0b)). Therefore, the jumps observed at the 
R3c —Pna2i (Oi) and the Pna2i (Oi) —Pca2i (O 2 ) 
structural transitions are directly linked to 9. 

In the Pnma phase, we observe two modes at 25 cm“^ 
and ~ 40 cm“^ (Fig0a)). The first one has quite large 
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spectral weight (10 times larger than all other magnon 
peaks), and does not seem to change frequency with in¬ 
creasing pressure. The second peak is quite similar to 
the magnon peaks in other phases. It may be that these 
first and second peaks are related to the so-called XS'*' 
and RS"*" antipolar modes known to occur in the non¬ 
polar Pnma phase Ea EO]. Future work is needed to 
verify such hypothesis because the model of Eqs. 0 -(i 
is presently developed to understand results of the polar 
R3c, Oi, O 2 and O 3 phases. 

In summary, via the combination of original Raman 
scattering experiment and computational and theoretical 
techniques, we report the magnetic excitations of BiFeOa 
as a function of pressure up to 12 GPa, showing for the 
first time that the material undergoes a series of mag¬ 
netic phase transitions linked to structural changes. As 
pressure increases above 3.5 GPa , the non-collinear cy¬ 
cloidal magnetism tranforms into a canted homogeneous 
antiferromagnet. At the ensuing crystal phase transitions 
at 5.5, 7.75 and 11 GPa, the two remaining spin excita¬ 
tions show jump discontinuities. Effective Hamiltonian 
approach provides information on the electrical polar¬ 
ization and structural changes of the oxygen octahedra 
through the successive structural phases. The extracted 
parameters are then used in a Ginzburg-Landau model to 
reproduce the evolution with pressure of the spin waves 
excitations observed at low energy in all the structural 
phases. We demonstrate that the structural phases and 
the magnetic anisotropy drive and control the spin ex¬ 
citations. Pressure reveals several hytherto unexplored 
regimes in the prototypical multiferroic BiFeOa and can 
help stabilize unstable structural distortions leading to 
promising novel metastable phases. The control of the 
crystallographic lattice by ultrafast optical excitations 
may result in high-speed magnonic thanks to the simul¬ 
taneous coherent driving of both lattice and magnetic 
excitations between different ferroelectric and magnetic 
phases. 
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